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Abstract 

Background: The host mounts an immune response to pathogens, but few data are currently available on the 
role of host genetics in variation in response to avian influenza (Al). The study presented here investigated the role 
of the host genetic background in response to in vivo infection with Al virus (A IV). 

Methods: Experimental lines of chicken and commercial crosses were experimentally infected intratracheal^ with 
10 3 EID 50 /bird of A/Chicken/ltaly/1 3474/99 H7N1 highly pathogenic avian influenza virus (HPAIV). Chickens were 
genotyped for the Mx polymorphism causing the S631N mutation, and for the Major Histocompatibility Complex 
(MHC). Whole-genome genotyping was carried out using 60 k Single Nucleotide Polymorphism (SNP) array 
developed by the poultry Genome-Wide Marker-Assisted Selection Consortium (GWMASC). 

Results: Variability in response of different chicken lines to the HPAIV infections and some degree of resistance to 
Al were observed: a statistically significant effect of chicken line on the response to infection was found. There was 
no association between survival in healthy conditions and polymorphisms at the Mx gene and the MHC-£ region. 
The analysis based on the 60 k SNPs provided a good clustering of the chicken lines, but no specific genetic 
cluster associated with response to AIV was identified. 

Conclusions: Neither the genotype at the Mx gene or MHC-£ locus, nor for SNP spanning the whole-genome 
identified loci involved in variations to response to AIV infection. These results point towards the possibility that 
either the genetic factors affecting the response of chickens to the H7N1 HPAIV are weak, or relevant alleles were 
not segregating in the studied populations. 



Background 

Avian influenza (Al) outbreaks are of concern as viral 
strains crossing the species barrier and infecting humans 
have been reported [e.g. 1]. The identification of genetic 
factors affecting resistance/susceptibility to pathogens 
would improve selection schemes for the development 
of resistant flocks of chicken. 

Host response to pathogens relies on innate and adap- 
tive immune responses. Mx proteins are antiviral mole- 
cules, which are part of the innate immune response 
[rev. 2], and genetic variations in the protein affect the 
antiviral activity. In chicken, Mx antiviral activity against 
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vesicular stomatitis virus and H5N1 virus has been 
reported to be dependent on the presence of asparagine 
at position 631 of the molecule [3]. 

Products of the Major Histocompatibility Complex 
(MHC) play a pivotal role in both innate and adaptive 
immune response, and chicken MHC haplotypes have 
been reported to be associated with either resistance or 
susceptibility to infectious diseases [rev. 4]. In order to 
identify additional loci that are involved in variations in 
response to infection, the chicken genome sequence [5] 
and Single Nucleotide Polymorphisms (SNPs) [6] infor- 
mation available make possible large scale and compre- 
hensive association studies to dissect genetic 
mechanisms behind the response of chickens to 
pathogens. 

Ltd. This is an open access article distributed under the terms of the Creative Commons 
s.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in 
properly cited. 
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In the present work, the effect of the genetic back- 
ground of a range of chicken lines on their response to 
in vivo AI virus (AIV) infection was investigated. 

Methods 

Experimental infection 

In addition to those previously described [7], six addi- 
tional broiler commercial lines, here referred to as lines 
F, G, H, I, L and M, were challenged with the A/ 
chicken/Italy/13474/99 H7N1 highly pathogenic avian 
influenza virus (HP AIV). Line F was a repetition of the 
experimental infection of line B, but with five more 
birds [7]. The experimental infection was carried out on 
20 9-week old birds for line F, and on groups of five 
adult birds each for lines G through M. The experimen- 
tal protocol was as described [7]. All the experiments 
were conducted complying to animal care committee at 
the institute where this study was carried out. 

Genotyping 

Genomic DNA was extracted from whole blood samples 
collected before the experimental infection. The 45 
birds constituting lines F to M were analysed for the 
Mx gene polymorphism causing the mutation at amino 
acid 631 as described [7]. 

The total of 130 chickens (including the 85 birds pre- 
viously described [7], and the 45 birds from lines F to 
M) were also genotyped for their variations at five 
microsatellite markers located within the MHC-£ 
region: LEI0258 [8], MCW0371 [9,10], MCW0312 [11], 
and two further microsatellite loci, MHC-D (F: 5'- 
CTGTTGGCGTTACAGAGCT-3'; R: 5'-TTCACCCAG- 
CAGCCTCTATC-3') and MHC-T (F: 5'-ATGGTGGC- 
CAAGTAAACTGGAG-3'; R: 5'- 

GGATCTGACAGCTGAGCGAGGT-3'). Each micosa- 
tellite marker was independently amplified. PCRs were 
carried out in a total volume of 12 with 30 ng geno- 
mic DNA, 1.2 |il lOx buffer, 1.5 mM MgCl 2 , 0.2 mM 
each dNTP, 3 pmol each primer, 0.5 U Taq polymerase. 
The PCR protocol was 94°C for 3 min, 30 cycles of 94°C 
for 30 s, annealing temperatures 56°C (MCW0312); 60°C 
(MCW0371 and MHC-D); 62°C (MHC-T); 65°C 
(LEI0258), for 30 s, 72°C for 30 s, and a final extension 
step at 72°C for 7 min. The forward primer of each 
microsatellite marker was 5'-labelled with either FAM or 
HEX. PCR products were analysed on ABI PRISM 3730 
DNA Analyzer (Applied Biosystems). Analysis of elec- 
tropherograms and allele scoring was performed with 
the GeneMapper software (Applied Biosystems). Haplo- 
types were reconstructed using the Arlequin software, 
ver 3.1 [12]. Statistical analysis was carried out using 
four inter-dependent variables (clinical score, time 
course of disease, chicken line, and MHC genotype) as 



previously described [7], and separately for the two vari- 
ables chicken lines and MHC genotypes. 

The 130 chicken samples were genotyped using the 60 
k SNP array developed by the poultry Genome-Wide 
Marker-Assisted Selection Consortium (GWMASC). 
Genotyping was carried out using the Illumina technol- 
ogy at DNA LandMarks Inc., Canada [13]. The genetic 
structure analysis was based on kinship of 17 chicken 
commercial lines. A total of 1,024 chickens were geno- 
typed. Fourteen samples randomly chosen were geno- 
typed twice as positive controls. GenABEL [14] was 
used to analyze the 60 k genotyping data within the R 
statistical environment [15]. The quality check on the 
dataset was carried out with the following thresholds: 
0.95 for marker call rate, 0.95 for sample call rate, 
0.0001 for minor allele frequency, as implemented in 
the "check.marker" function of GenABEL. Samples shar- 
ing genotype identity higher than 0.95 with another 
sample were excluded from the analysis. No markers 
were excluded following the test for Hardy- Weinberg 
equilibrium. Classical Multidimensional Scaling based 
on kinship analysis of the clean dataset was used to 
define samples in two-dimensions, and clusters were 
defined "a posteriori" based on kinship analysis. 

Results and discussion 

Clinical observations 

All birds from lines H, I, L, and M died within 6, 9, 6, 
and 10 days after the infections (PI), respectively (Fig. 
1). Only 2 birds of line G survived the infection (the 
remaining 3 birds of line G died within day 5 PI). Line 
F was a repetition of line B [7]. All birds from the for- 
mer group died within 5 days following infection. 

Genetic analyses 

The six lines (F to M) were analysed for the G/A poly- 
morphism at position 2,032 of the chicken Mx cDNA 
and were segregating for the two alleles, as previously 
observed for the broiler commercial lines B and C [7]. 
The A allele, coding for asparagine at position 631 of 
the protein putatively associated with higher antiviral 
activity [3], had frequencies ranging between 10% (lines 
I and M) and 90% (line G). The statistical analyses car- 
ried out on all 130 birds (including all the birds from 
lines A to E [7]) did not find an association between the 
Mx genotype and the response of the birds to AIV 
infection as previously found by Sironi and colleagues 
[7]. However, a significant effect of the line on clinical 
status after the viral challenge was found [7], thus sug- 
gesting a genetic basis to variability in response to the 
pathogen. 

Twenty-eight MHC haplotypes reconstructed using 
the genotyping data across the five microsatellites were 
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Figure 1 Survival plot. Survival curves for 6 chicken lines following experimental challenge (day "0") with HPAI H7N1 influenza virus. 



identified among the birds used in the study, constitut- 
ing a total of 60 genotypes. The genotype frequency ran- 
ged from 0.77% to 11.5%, and 37 out of the 60 
genotypes were unique (i.e. present only once in the 
population analysed). Higher homozygosity at the 
MHC-B region was observed in lines E and B compared 
to the other lines, confirming that these lines contained 
greater genetic diversity as observed for genotyping 
results obtained using microsatellite markers scattered 
through the genome (data not shown). Twelve of the 60 
MHC genotypes were found among the birds that sur- 
vived the infection: four birds carried unique MHC gen- 
otypes; whereas the remaining 18 surviving birds carried 
genotypes also present in birds that died as a result of 
the AIV infection. The statistical analysis, in which the 
effect of the two variables "line" and "MHC genotype" 
were estimated independently, showed that the MHC 
genotype did not have significant effect on the survival 
in healthy clinical status (p = 0.3574). Further associa- 
tion studies based on chicken groups characterised by 
specific serotypes are fundamental to understand if 
MHC plays a significant role in response and outcome 
of infection with AIV. 

Kinship analysis of genotype data from the 60 k SNP 
set clustered the majority of birds into the 11 lines stu- 
died. All birds from lines A, E, I, and M clustered into 
separate unique clusters, whereas line C birds split into 
2 clusters. However, the close proximity of the clusters 
suggested highly similar genetic background among 
these chickens. The same was observed for lines G and 
H, each of which split into two very closely related clus- 
ters. Additionally the clusters from these lines were in 
proximity indicating that they had similar genetic 



background. A group of birds from line C overlapped 
with the cluster for line D. Finally, lines B, F, and L 
clustered together. This was expected for lines B and F: 
the latter to be an experimental repetition of the former, 
even though chickens from lines B and F showed differ- 
ent responses to AIV infection. In order to further 
investigate a possible genetic distinction between lines B 
and F, an additional analysis was carried out with GenA- 
BEL using exclusively the samples from these lines. 
Using all permutations of possible cluster definitions 
(pre-setting 2, 3, 4, and 5 as potential number of clus- 
ters) individuals from the two lines were split over all 
available clusters, indicating that lines B and F were 
indeed genetically closely related. 

Conclusions 

This study showed that genetic line has a statistically 
significant effect on the response of chickens to the 
H7N1 HPAIV. Nevertheless, no significant association 
between variability in response and polymorphisms in 
specific candidate loci (Mx and MHC-.S) or genome- 
wide (60 k SNP array) was identified. 

As the Mx gene was only genotyped for the 2,032 
polymorphism, it cannot be ruled out another poly- 
morphism, within the Mx gene, may influence the in 
vivo antiviral activity. It is also possible that alleles of 
other genes that were not segregating in the lines exam- 
ined in the present study confer resistance to H7N1 
HPAIV. 
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